AD-A098  476  VIRGINIA  POLYTECHNIC  InST  AND  STATE  UNIV  BLACKSRURG  D— ETC  P/G  12/1 
ADAPTIVE  ARMA  SPECTRAL  ESTIMATION* (U) 

1981  J  A  CADZOW*  K  OGINO  N00014-80-C-0303 

UNCLASSIFIED  4  NL 


AD  A098476 


sscurity  classification  of  this  rags  r»in  Qtm  &nn< 

I  REPORT  DOCUMENTATION  PAGE 


PQftT  NUMICR 


T.  JSPQIIT  N 

m  \  /■ 


Tl.  OO VT  ACCCUtOM  HO. 


«.  TITLC  (m*4  Su*lltf) 

\  Adaptive  ARMA  Spectral 


IJ  Estimation^ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
i.  hscificnt’s  catalog  humScr 


S.  TYR*  or  RSRORT  i  RCftlOO  COVCRCD 

Manuscript 

I  rcnrorming  ora.  report  mumgcr 


[7.  AyiMoiw 

James  A.  jc adzow 
Koj ijoglno 


1.  contract  on  on *nt  NUMee*r«j 


*•  RCRRORMING  ORGANIZATION  NAMC  ANO  AOORef^-..  - 

Virginia  Tech  -  —  1  Jl'X  l 

Department  of  Electrical  Engineering  { 

Blacksburg,  VA  24061 


11.  CONTROLLING  ORFICC  NAMC  ANO  AOORCSS 

Office  of  Naval  Research  (Code  436 )( //, 
Statistics  and  Probability  Program  V_  ■ 
Arlington.  VA  22217 


4.  MONITORING  AGCNCY  NAMC  A  AOORCSV"  dllltmtt  I 


12.  ACPOMT  OATC 

13.  NUMtCRO®  PACES 

IS.  SECURITY  CLASS.  (•!  t fit* 

Unclassified 

rmpowt) 

ts«.  OCCLASSIPICATION/  DOWNGRADING 
SCHEDULE 

17.  OlSTRiauTlON  JT  ATSMCNT  (9 1  1*9  9»9tr9 cl  m 99994  in  *(•«*  70,  II  dlll9t9lH  from  lifMJ 


I*.  iurrlcmcntary  notes 


WT®  BOCWSJRST  IB  BEST  QUALITY 

***  OQPg  njRSIISRBD  TO  DCC  A~ 

UGBIFICAJff  NUMBER  07  PAGJL3  UHICBi  DO  SOS 


ft  Key  WO  AOS  tfn  </  •!•«*  wwiH&MJ 

ARMA  Model,  Adaptive  Pararaeteric  Method,  Yule-Walker  Equations 


ip.  RACT  rComlms*  «i  KM  1/  rbcmmvb  m*W  Ry  Mm®  immms) 

JA  novel  adaptive  method  for  efficiently  obtaining  an  ARMA  model  spectral  esti  late 
of  a  wide-sense  stationary  time  series  is  presented.  It  is  adaptive  in  the  s  inse 
that  as  a  new  element  of  the  time  series  is  observed,  the  coefficients  of  a  (>,p)ch 
order  ASHA  model  may  be  algorithmically  updated.  This  algorithm’s  computatio  lal 
complexity  (i.e.  the  number  of  multiplications  and  additions  required)  is  of  :he 
order  p  log(p)  for  a  particular  version  of  the  method.  Moreover,  the  spectra  . 
estimation  performance  of  this  new  method  is  found  typically  to  be  far  super! >r 
to  such  contemporary  approaches_as_the_Box-Jenkins1_maximum_entrogy ,  _  ^ 

00  ,  1473  coi  tion  or  ,  NOV  .»  If  ogk>lct«  Unclassified  ^  X  J-0  6>  ^ 

i  N  3102'..  •  2 1 J ■  3  - -- 


Unclassified  J-  u>  <j 

siclnity  :l AftirtCATicN  o *  *»ii  rags  *n«r  i»«mi 


/  ,  i  f  i 


l 


SECURITY  CLASSIFICATION  OF  THIS  RACE  (•*>«■  Om»  tfXf*0 


U20.  (ConC'd) 

jWidrow's  LMS  methods.  This  performance  in  conjunction  with  its  computational 
efficiency  mark  this  algorithm  as  being  a  primary  spectral  estimation  tool. 

'  4^ 


S/N  3 1 02-  Lf-0w.««01 


"Accession  For. 

"UTIS  GRAM 

DTIC  tab 

Unannouncort 

justi fioatior 

•Qy__ _ 

pistritution/ 

fAvaiieMUtY  *£• 

1  .  Avail  anc/o. 

'  Saocial 

■plot 


SECURITY  CLASSIFICATION  OF  Tt»u  RAC! trnm  Omm  ImmR 


v* 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


ADAPTIVE  ARMA  SPECTRA!  ESTIMATION 


James  A.  Cadzow  &  Koji  Qgino 
Department  of  Electrical  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 


ABSTRACT 

A  novel  adapcive  method  for  efficiently  obtain¬ 
ing  an  ARMA  model  spectral  estimate  of  a  wlde-sense 
stationary  time  series  Is  presented.  It  is  adapc¬ 
ive  in  the  sense  chat  as  a  new  element  of  the  time 
series  is  observed,  the  coefficients  of  a  <p,p)ch 
order  ARMA  model  may  be  algorithmically  updated. 

This  algorithm’s  computational  complexity  (i.e., 
the  number  of  mulclplications  and  additions 
required)  is  of  the  order  p  log(p)  for  a  particular 
version  of  the  method.  Moreover,  che  spectral 
escimacion  performance  of  this  new  method  is  found 
typically  to  be  far  superior  to  such  contemporary 
approaches  as  che  Box-Jenkins,  maximum  entropy, 
and,  Widrow's  LMS  methods.  This  performance  in 
conjunction  with  its  computational  efficiency  mark 
chls  algorithm  as  being  a  primary  spectral 
estimation  tool. 


I.  INTRODUCTION 


In  various  signal  processing  applications,  it 
is  necessary  to  escimate  the  power  spectral  density 
of  a  vide-sense  stationary  time  series  {x„}.  Since 
only  a  finite  sec  of  time  series'  observations  are 
typically  available  for  this  cask,  one  almosc 
ilvays  invokes  a  finite  parameter  model  for  the 
soeocral  density  escimace.  Vlchouc  doubc,  che 
rational  spectral  density  function  as  specified  by 
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constitutes  the  most  widely  used  of  such  models. 
This  particular  model  is  generally  referred  to  ae 
being  an  autoregressive-moving  average  (ARMA)  model 
of  'rcer  <p,q). 

The  preaominanc  effort  In  rational  spectral 
estimation  has  been  directed  towards  the  more 
specialized  aucoregressivc  (AR)  model  for  which 
1  •  !.  and.  Cne  moving- average  <MA)  model  in  which 
1  •  l.  As  examples,  the  maximum  entropy,  one-step 
■radiator,  ano  lutoregressive  methods  have  been 
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developed  for  efficiently  estimating  the  AR  model's 
Afe  coefficients.  Similarly,  the  perlodogram  and 
its  variants  have  been  found  to  yield  effective 
MA  modeling  procedures.  The  Interested  reader  will 
find  excellent  treatments  of  these  and  ocher 
rational  spectral  estimation  methods  in  Haykin  [1] 
and  Childers  [2]. 

Recently,  attention  haa  been  focused  on 
developing  so-called  super-efficient  algorithms  for 
estimating  the  AR  model's  a^  coefficients  whereby 
on  the  order  of  p  log(p)  computations  are  required 
for  this  task.  These  algorithms  are  typically 
predicated  on  the  divide  and  conquer  approach 
(e.g.,  see  refs.  [  3 ]&[  4 ]).  These  super  algori¬ 
thms  offar  the  potential  of  providing  a  significant 
computational  advancaga  in  ganarating  AR  spectral 
estimates  when  compared  to  ocher  contemporary  AR 
procadurae.  Unfortunately,  implementation  of  these 
super  algorithms  is  relaciveiy  complex  and  a  rather 
large  value  for  the  AR  order  parameter  p  le 
required  before  che  computational  complexity 
p  log(p)  is  approached.  It  is  felt  Chet  future 
developments  will  alleviate  these  difficulties. 

Despite  the  concentrated  interest  given  to 
AR  spectral  escimacion.  it  is  widely  recognized 
chat  an  ARMA  spectral  model  is  generally  the  most 
effective  rational  model  from  a  parameter  parsi¬ 
mony  viewpoint.  In  recognition  of  this  fact,  a 
variety  of  procedures  have  been  developed  for 
generating  ARMA  models.  These  include  che  whiten¬ 
ing  filter  approach  which  is  typically  iterative 
in  nature,  generally  slow  in  convergence,  and. 
usually  requires  an  excessively  large  number  of 
time  series'  observations  to  be  effective  (e.g., 
sea  refs.  [3]  6  [6]).  More  desirable  closed  fora 
procedures  which  overcome  chase  deficiencies  have 
been  offered.  These  Include  che  so-called 
3ox-Jenklns  method  and  ics  variants  [  7  1  --  [91. 
and,  more  recendy,  Cadzow  has  developed  a  "high 
performance"  mechod  [101  4  [Hi-  Although  this 
laccer  mechod  haa  provided  excellent  speccral 
estimation  performance  when  compared  to  che 
maximum  entropy  and  3ox-Jenkins  methods,  its 
:omputationai  efficiency  is  somevhac  inferior. 

..'e  shall  herein  presenc  i  novel  ilgebraic 
joprcach  for  generating  in  ARMA  modal  spectral 
estimate.  It  offers  the  dual  advancaga  jf  having 
a  suaer  ilgorlthm's  computational  efficient-'  wniie 
ac  the  same  time  maintaining  a  ssactral  estimation 


capability  similar  to  che  abova  mentioned  high 
performance  method.  These  characteristics  mark 
chis  algorithm  as  being  a  primary  spectral  estima¬ 
tion  tool. 

In  this  paper,  we  shall  consider  exclusively 
che  task  or  estimating  che  ARMA  model's  auto¬ 
regressive  coefficients.  This  estimation  is,  to 
a  large  extenc,  motivated  by  che  well-known  Yule- 
Walker  equations  which  will  be  briefly  reviewed 
in  the  next  section.  Once  these  aucoregressive 
coefficient  estimates  have  been  obtained,  a 
variety  of  procedures  exist  for  estimating  the 
ARMA  model's  moving-average  coefficients.  The 
interested  reader  may  consult  references  [10]  6 
[11]  for  a  description  of  such  procedures. 


II.  FUNDAMENTAL  CONCEPTS 


It  is  readily  shown  chat  che  random  time 
series  wich  spectral  densicy  (1)  can  be  modeled  as 
being  che  response  of  the  causal  ARMA  system 
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(1) 


to  che  zero  mean  whice  noise  excitation  (en;  whose 
individual  terms  have  variance  :2,  Although  the 
more  general  case  may  be  straightforwardly  created, 
we  have  here  restricted  q  ■  p  for  purposes  of 
ease  of  understanding.  The  autocorrelation 
characterization  of  chis  ARMA  system  is  readily 
achieved  by  first  multiplying  each  side  of  equaeion 
(2)  by  the  entity  xjj..,,,  and  then  caking  che  expect¬ 
ed  value.  This  results  in  che  well  known  Yule- 
Walker  equations  as  given  by 


corresponding  to  che  enlarged  n+1  data  sec 
(i.e.,  x1>  x,,  .  .  .  ,  xQ+l)  . 

III.  ARMA  ESTIMATION:  DIRECT  APPROACH 

In  chis  section,  a  procedure  for  estimating 
che  ARMA  models'  aucoregressive  coefficients  shall 
be  given.  To  begin  chis  development,  one  first 
evaluates  Che  model  equation  (2)  over  the  sec 
p  +  2  1  k  ^  a  to  obtain  the  following  time  series' 
relationships 
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It  will  be  convenient  to  represent  this  relation¬ 
ship  in  che  more  compact  vector  format 

x  X  a  •  £  b  :6b) 


P 

r  (m)  +•  7  a.  r  (m-k)  *  0  for  m  >  p  (3) 

X  ,  - .  (V.  X 

In  chis  equaeion,  che  symbol  rx(m)  denotes  the 
time  series'  autocorrelation  sequence 


where  *  and  E  denote  che  operations  of  complex 
conjugation  and  expected  value,  respectively. 

In  what  is  to  follow,  the  Yule-Walker  equations 
1 will  serve  as  a  motivating  influence  In  evolv¬ 
ing  a  method  for  estimating  the  aucoregressive 
-.’efficients  of  che  ARJIA  model  (2).  These  auto¬ 
regressive  coefficient  estimates  are  to  be  based 
totally  on  che  following  contiguous  set  of  n  time 
series  observations 


The  method  to  be  described  will  be  adapeive  In 
nature.  Samel'.",  as  the  new  time  series  element 
■.n-:  becomes  avatlaoie.  It  is  possible  to 
efficiently  uoaate  che  opclmal  aucoregressive  co¬ 
efficients  generated  from  the  n  data  set  (3)  to 
brain  the  ’petmai  aucoregressive  coef f iclents 


where  x,  a,  and  ]>  are  n-p-1,  p.  and,  p+1  column 
vectors  respectively,  while  X  and£  are  '.n-p-1)  ■? 
and  (n-p-1)  »  (p+1)  matrices,  respectively.  The 
encries  of  these  vectors  and  matrices  are  obtained 
by  directly  comparing  relationships  (6a)  and  (6'o). 

We  now  wish  to  use  che  Yule-Walker  equations 
(3)  In  conjunction  wich  relationship  (6)  to 
estimate  the  aucoregressive  coefficient  vector  a. 
This  objective  is  readily  achieved  by  first  intro¬ 
ducing  the  following  (n-p-1)  *  c  lower  triangular 


type  matrix 
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- 

- 

the  selection  or  che  integer  c. 

whicn  specifies 

the  number  of  columns  of  matrix  Y  ,  is  critical. 
A  discussion  of  how  one  goes  about  making  this 
selection  will  be  shortly  given. 

Upon  left  multiplying  each  sice  of  reiacion- 


i 


ship  (6b)  by  eh*  complex  conjugate  transpose  of 
matrix  (7)  as  denoted  by  Y*  ,  chars  results 

Yrx  +  Y~X  a  -  Y*£  b  (8) 

this  syscaa  or  equations  is  raadlly  found  co  con- 
sclcuca  a  statistical  approximation  to  tha  first 
t  Yula-Walkar  aquations  (i.a.,  axprasslon  (3)  for 
p  <  a  ^  p+t).  This  is  raadlly  varlfiad  by  caking 
che  expected  value  of  expression  (8)  which  results 
in  che  following  sec  of  aquations 

P  -  ] 

(n-m)  r  (m) +•  J  a^r^m-k)  «0  for  p<a<p+c 

L  k*1  J  "  (9) 

It  is  co  be  noced  that  the  right  side  zero  Cera 
arises  due  to  che  face  chac  E(Y+£;  *  3  (i.e.,  the 
null  matrix).  This  la  a  direct  consequence  of  the 
A3MA  models'  causality  and  the  whiteness  of  the 
excitation  time  series  which  causes 

E{x  *  s,  }  -  0  for  k  >  n. 
n  k 

With  che  above  observations  in  mind,  a  logical 
selection  procedure  for  the  autoregressive  co¬ 
efficient  vector  is  suggested.  Namely,  che  auto¬ 
regressive  coefficient  vector  a  is  selected 
so  as  co  cause  che  left  side  of  relationship  (8) 
to  be  as  close  as  possible  co  che  zero  vector 
(i.e.,  the  expected  value  of  the  right  side  vector 
Y~r  b  ) .  This  results  in  an  approximation  to  the 
Yule-Walker  equations  which  Is  "most"  consistent 
with  che  time  series'  observations  (5).  A  parti¬ 
cularly  convenient  measure  of  che  closeness  of 
YV  b  to  the  zero  vector  is  given  by  the  quadratic 
functional 

f(a)  -  (Y~x  +■  Y*X  a  )1’W(Y',’x  +  Y*Xa)  (10) 

in  which  W  is  a  c  »  t  symmetric  positive  Semi¬ 
te:  inice  matrix  chat  is  usually  selected  so  as  co 
weight  differently  various  elements  of  the  error 
vector  Y  x  +  Y*Xa.  Ic  is  readllv  shown  that  the 
autoregressive  coefficient  vector  which  minimizes 
this  quadratic  functional  satisfies  che  following 
nsiscenc  system  of  p  linear  equations  in  the 
?  autoregressive  coefficient  unknowns 

X'YWY'Xa’  «  -  x'VVY'x  (11) 

Ic  is  possible  co  use  projection  theory  concepts 
to  achieve  a  oompucac lonally  efficient  method  for 
'bcuintr.g  the  optimum  vector  a*  when  che  weight¬ 
ing  matrix  N  is  exponentially  diagonal.  A  paper 
now  in  preparation  will  detail  this  solution  pro- 
--'jure. 

A  few  words  ire  now  appropriate  concerning 
cna  selection  of  the  integer  t  which  in  part 
-xuracterlzes  matrix  Y  as  given  by  expression 

If  t  is  sec  equal  to  p  ,  it  is  seen  that 
relationship  .  7)  constitutes  a  statistical 
estimate  if  the  first  p  Yule-Walker  equations, 
is  suen,  reiaevonshio  ill)  bears  a  resemblance  to 


the  Box-  Jenkins  sechod  of  autoregressive  co¬ 
efficient  eetlaetlon  (  7  ) .  Upon  cloeer  examination, 
however,  it  It  found  that  thnee  two  approaches  are 
quite  different.  As  a  setter  of  fact.  It  has  been 
empirically  found  that  che  spectral  estimation 
performance  which  results  from  utilization  of 
raletlonahip  (11)  with  c  •  p  Is  distinctly  better 
than  that  obtained  with  the  Sox- Jenkins  aaehod.  [11]. 

If  che  Integer  t  Is  taken  co  be  larger  chan 
p,  chan  eore  than  che  alniaal  number  (i.e.,  p)  of 
Yule-Ualker  aquation  approxiaacions  are  generated. 
With  this  larger  base  of  Yule-Uelker  equation 
approximations  (1. a.,  greater  than  p) ,  Ic  night  be 
conjectured  chat  an  iaprovaaent  In  spectral  esti¬ 
mation  performance  would  result.  This  anticipated 
Iaprovaaent  In  perforaence  has  been  In  feet 
empirically  demonstrated  on  numerous  examples 
craatad  co  data. 

It  is  raadlly  shown  that  che  procedure  here 
presented  for  selecting  che  autoregressive  co¬ 
efficients  is  tquivalent  to  Cadzow's  high  perform¬ 
ance  AftMA  tpeccrsl  estimation  method  [10]&[11J.  As 
such,  a  large  bate  of  empirical  evidence  gathered 
in  using  this  latter  procedure  suggests  chat  che 
spectral  estimation  performance  of  this  paper's 
proesdura  is  clearly  superior  to  chac  achieved  by 
such  commonly  used  contemporary  procedures  as  the 
Box- Jenkins  and  maximum  entropy  methods.  The 
advantage  accrued  in  algebraically  formulating  che 
spectral  estimation  problem  via  relationship  (8) 

(as  comparad  to  che  equivalent  high  performance 
method)  resides  in  the  ability  co  directly  use  . 
sophisticated  least  mean  square  concepts. 

Using  these  concepts,  it  is  possible  co  evolve  a 
computationally  efficient  adaptive  spectral  esti¬ 
mation  algorithm  whose  computational  complexity  is 
on  the  order  of  p^  .  More  importantly,  by  an 
appropriate  modification  of  the  vector  x  and 
matrices  X  and  Y  used  in  relationship  (8),  it 
is  posslbls  co  obtain  a  truly  super  computational 
algorithmic  procedure  which  has  a  computational 
complexity  of  che  order  p  log(p).  It  is  this 
capability  which  distinguishes  the  herein  presented 
method  from  che  high  performance  method. 


IV.  AKMA  ESTIMATION :  MODIFIED  APPROACH 

It  is  possible  co  realize  a  significant  compu¬ 
tational  Improvement  In  che  proposed  spectral 
astimacion  algorithm  by  rastrlctlng  c  to  be  ? and 
by  appropriately  modifying  Che  vector  and  matrix 
entrees  in  expression  (8).  Although  e  varisty  of 
such  modifications  art  faaslbia.  we  will  here 
reecrlcc  our  intereets  co  two  such  possibilities 
calisd  ehs  "premodif lcaclon”  and  ''postmodificacion'' 
approaches.  Additional  modifications  are  being 
examined  and  will  be  shortly  reported  upon. 

(a)  Premodification  Method 

In  the  premodification  mechoq.  the  x  vectqr 
and  V  matrix  as  given  in  exprtsslons  <6>  and  \"i 
remain  the  same  vnlle  the  integer  t  is  fixed  at 
?.  The  matrix  X.  however,  is  modifisd  to  che 


I 


following  low«r  triangular  format 


•  0  I 


If  this  matrix  la  substituted  into  relationship 
(3) ,  an  alternate  method  for  estimating  the  auto¬ 
regressive  coefficients  Is  ac  hand.  Due  to  the 
lower  triangular  structure  of  boch  matrices  X  and 
Y,  however.  It  Is  possible  to  Implement  an  adapt¬ 
ive  algorithm  for  obtaining  such  escimatas  with  a 
computational  complexity  of  order  p.  This  remark¬ 
able  Improvement  in  computational  efficiency  Is 
found  to  be  dependent  on  a  lattice  structure  Imple¬ 
mentation  of  the  ARMA  model  (e.g.,  see  ref.  [13]). 

A  measure  of  the  spectral  performance  of  this 
preaodif ication  method  may  be  obtained  by  taking 
the  expected  value  of  expression  (8)  with  the 
matrix  substitution.  This  Is  readily  found  to 
yield 


(n-mi  _  a.^r  ^m-k)  *■  _  (n-p-k-l)a.  r  (m-k)  « 0  (13) 

k«0  k*m-p  ,r' 

p  <  m  Op 

where  tne  aq  coefficient  is  set  equal  to  one. 
Although  these  relationships  don't  precisely  satis¬ 
fy  the  Yule-Walker  equations  as  in  the  direct 
ibproach.  it  Is  noted  mat  for  n  >>  p,  an  excell¬ 
ent  lpproxiavation  is  In  fact  realized.  With  this  ' 
m  mine,  it  is  reasonable  to  anticipate  that  this 
iremodi: ication  method  will  have  a  apectral  escl- 
nocion  oerformance  approaching  thac  of  the  afore¬ 
mentioned  hlah  Performance  (direct)  method.  Em- 
"'‘.ricil  eviaence  cathered  to  date  re-enforces  this 


3  stmoalf itjtion  Method 

In  tne  oostmodif icat ion  memod,  the  altera- 
:  i. ns  to  oe  made  are  given  by 


<n-l  "-n  0 . 0 


where  the  prime  denotes  the  transpose  operator. 
Using  these  entries  in  expression  (3) ,  it  is  possi¬ 
ble  to  evolve  an  adaptive  algorithmic  solution  pro¬ 
cedure  for  the  optimum  autoregressive  coefficients. 
This  entails  utilization  of  the  doubling  algorithm 
concept  and  results  In  an  adaptive  algorithm  whose 
computational  complexity  is  on  the  order  of 
P  log(p)  [13]. 

To  gauge  the  effectiveness  of  this  postmodifi¬ 
cation  method  as  embodied  In  expression  (8),  the 
expected  value  of  this  expression  with  entries  (14) 
is  next  taken  and  results  In 

P 

(n-m)rx(m)  ♦  ^  (a-»Hc)a.  r  (m-k)  •  0  p<a<_2p  (13) 

k*l 

As  In  the  premodification  method,  these  expected 
value  relationships  do  not  precisely  satisfy  the 
governing  Yule-Walker  equations.  On  the  other  hand, 
for  a  >>  p,  It  is  clear  chat  they  do  provide  a  very 
excellent  approximation  to  these  characteristic 
equations.  As  such,  it  is  not  surprising  chat  Che 
autoregressive  coefficient  estimates  provided  by 
the  poscmodlficaclon  method  as  represented  by 
expression  (3)  yield  a  very  satisfactory  spectral 
estimation  performance.  This  performance  has  been 
empirically  found  to  approach  chat  of  che  high 
performance  (direct)  method. 

In  order  to  cast  the  effectiveness  of  che 
herein  proposed  AKMA  spectral  estimation  approach, 
cha  classical  problem  of  resolving  cwo  sinusoids 
ambedded  In  additive  white  noise  was  considered. 

The  spectral  asclmaces  obtained  from  these  methods 
are  shown  In  Fig.  1  along  with  che  results  genera¬ 
ted  using  the  maximum  entropy  and  Box-Jenkins 
methods.  It  is  apparenc  thac  the  estimates  obtain¬ 
ed  using  this  paper's  approach  were  clearly  super¬ 
ior  for  this  cask. 

V.  CONCLUSION 

A  novel  approach  to  ASMA  model  spectral  esti¬ 
mation  haa  baan  prasanted.  This  esclcaclon 
approach  posaasaaa  tha  dual  atcrlbutaa  of  providing 
an  exeallant  spactral  estimation  performance,  and, 
of  having  an  exceptional  computational  efficiency. 
Its  excellent  spactral  estimation  performance  has 
bean  demons cratsd  on  numarous  examples  treaced  to 
decs  and  virtually  always  exceeds  chat  achieved  by 
such  conesmporary  procedures  as  tha  maximum  entropy 
and  Box-Jenkins  Methods. 


Tha  above  mentioned  super  computational  effic¬ 
iency  wee  achieved  by  appropriately  modifying  some 
vector  and  matrix  entrlee.  In  particular,  cwo  such 
modifications  have  been  herein  offered.  Further 
studies  are  now  being  conducted  relative  to  employ¬ 
ing  cha  basic  approach  hareln  caken  to  evolve  even 


batter  performance.  These  results  will  be  reported 
In  forthcoming  papers. 
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